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LONG-TERM  GOALS 

This  is  part  of  a  collaborative  program  Optimal  Deployment  of  Drifting  Acoustic  Sensors,  ODD  AS, 
which  involves  academic,  private  sector,  and  Navy  scientists.  The  long-term  goal  of  ODD  AS  is  to 
develop  strategies  for  the  deployment  of  drifting  sensors  that  maximize  the  amount  of  environmental 
information  collected  with  the  fewest  sensors.  The  primary  long-term  goal  of  our  component  is  to 
quantify  oceanic  submesoscale  stirring  and  advective  transport  processes  and  to  assess  their  roles  in 
the  deployment  strategies.  The  secondary  goal  is  to  develop  novel  and  quantitative  assessments  of  the 
ability  of  predictive  models  and  real-time  observations  to  determine  the  movement  of  sensors  and  the 
location  of  material  boundaries  or  advective  pathways  in  the  ocean. 

OBJECTIVES 

Two  objectives  were  pursued  during  EY06: 

•  Locate  hyperbolic  trajectories  and  calculate  associated  inflowing  and  outflowing  manifolds 
in  EAS16  hindcasts  for  the  region  around  Taiwan. 

•  Assess  the  ability  of  EAS16  to  predict  the  movement  of  drifting  sensors. 

APPROACH 

As  documented  in  publications  and  previous  progress  reports  we  have  developed  a  variety  of 
Lagrangian  methods  that  quantify  advective  processes  from  synoptic  velocity  archives.  These  methods 
are  used  in  this  project  to  configure  deployment  schemes  of  drifting  sensors  so  as  to  maximize  data 
return  from  drifting  sensor  arrays  with  a  minimum  number  of  units.  This  requires  meeting  three 
diverse,  and  occassionally  contradictory,  user  requirements.  One  is  to  maximize  the  spatial  coverage 
for  regions  of  interest.  Another  requirement  is  maintain  array  coherence  for  a  specified  time  period. 
The  third  requirement  is  to  sample  denied  areas.  To  meet  these  requirements  it  is  also  important  to 
decide  whether  to  deploy  all  sensors  at  one  time  or  to  use  a  time  staggered  strategy. 

The  key  of  our  approach  is  the  capability  to  locate  critical  regions  from  the  synoptic  velocity  archives. 
These  are  usually  found  in  regions  where  the  velocity  is  small.  However,  small  velocity  regions  in  the 
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ocean  have  vastly  different  characteristics.  Hyperbolic  regions,  typieally  found  in  the  nether  region 
between  oeean  eddies,  are  eharacterized  by  strong  shear.  Deployment  in  these  regions  is  appropriate 
for  broad  area  data  coverage  and  possibly  aceess  to  denied  regions.  However,  hyperbolic  regions 
should  be  avoided  when  the  primary  goal  is  array  coherence  sinee  these  regions  typieally  have 
extremely  strong  shear  that  produces  extensive  array  distortion.  See  Toner  and  Poje  (2004)  for 
additonal  discussion  of  the  importance  of  hyperbolie  regions  on  deployment  strategies. 

The  data  used  in  this  effort  comes  from  two  sources.  The  model  results  come  from  the  Navy  EAS16 
hindeasts  for  the  region  around  Taiwan.  Drifting  sensor  trajectories  eome  from  deployments  by  SIO 
and  NAVAIR  exercises.  For  the  first  projeet  objective,  EAS16  velocities  are  used  to  determine 
hyperbolie  regions  in  the  planned  deployment  areas  and  to  assess  the  evolution  of  hypothetieal 
deployment  arrays.  For  the  seeond  objective  we  compare  the  predieted  movements  of  individual 
drifting  sensors  from  EAS16  with  the  observed  trajectories.  Two  eriteria  are  available  for  quantitative 
assessment.  These  are  the  Eagrangian  error  metries  proposed  by  Toner  et  al.  (2001),  and  the  aeeuracy 
of  advective  pathways  predicted  by  inflowing  and  outflowing  manifolds  assoeiated  with  hyperbolic 
trajectories. 

The  researeh  on  this  effort  is  the  result  of  collaboration  with  A.  Poje  at  CETNY  Staten  Island,  R.  Allard 
and  S.  Riedlinger  at  NRL  Stennis,  the  drifter  group  at  SIO  led  by  P.  Niiler,  D.  Fenton  and 
D.  Hammond  at  NAVAIR,  and  R.  Heitsenrether  at  APE  JHU.  The  NRL  group  supplies  model 
hindeasts  from  EAS16,  NAVAIR  and  JHU  supply  data  from  fleet  exercises.  The  primary  responsibility 
of  the  group  at  the  University  of  Delaware,  working  with  A.  Poje,  is  the  Eagrangian  analysis. 

WORK  COMPLETED 

During  this  fiscal  year  we  concentrated  on  finding  ways  to  loeate  hyperbolic  regions  and  calculate 
manifolds  from  the  EAS16  model  velocity  fields.  The  problem  is  acute  since  the  model  fields  contain  a 
strong  tide  and  inertial  signal.  These  effeets  mask  the  separation  of  Eulerian  and  Eagrangian  time 
scales  neeessary  for  present  methods  to  work.  Haller  and  Poje  (2001)  and  Kirwan  (2006)  provide 
detailed  teehnical  discussions  of  this  point.  We  had  some  sueeess  in  loeating  these  regions  and 
growing  the  manifolds  by  working  with  daily  averaged  fields.  The  manifolds  were  then  superposed  on 
the  hourly  fields.  Comparisons  with  model  trajectories  calculated  from  the  hourly  fields  were  then 
made.  Preliminary  results  are  summarized  in  the  next  section. 

Hurricane  Katrina  slowed  aeeess  to  EAS16  model  velocity  fields. 

RESULTS 

Hyperbolie  trajeetories  are  located  at  the  intersection  of  two  material  manifolds.  These  manifolds 
appear  as  eurves  on  planar  oeean  surfaees,  but  in  reality  they  are  two  dimensional  fronts  that  extend 
well  into  the  water  column.  They  exist  for  finite  periods  of  times,  often  of  the  order  of  a  week  or 
longer.  During  this  period  they  may  extend  aeross  several  hundred  kilometers  of  oeean.  Sinee  they  are 
impervious  to  adveetion,  they  provide  a  qualitative  template  for  parsing  the  origins  and  fates  of 
partieles.  Fluid  near  the  inflowing  or  “stable”  manifold  is  drawn  toward  the  hyperbolic  trajectory  while 
fluid  near  the  outflowing  or  “unstable”  manifold  is  repelled.  Partieles  near  each  other,  but  on  opposite 
sides  of  an  inflowing  manifold,  will  diverge  exponentially  in  time  as  they  approaeh  the  hyperbolic 
trajectory.  As  noted  by  Haller  and  Poje  (2001)  this  produees  a  separation  of  Eagrangian  and  Eulerian 


timescales  in  hyperbolic  regions.  We,  and  others,  have  used  this  timescale  separation  effectively  in 
locating  hyperbolic  trajectories. 

Unfortunately,  strong  tides  and  inertial  flows  affect  both  the  hyperbolic  trajectory  and  the  nearby 
particles,  thus  camouflaging  the  time  scale  separation.  This  has  prevented  us  from  finding  hyperbolic 
trajectories  in  the  EAS16  model  field  in  the  East  Asia  Sea.  During  this  fiscal  year  we  addressed  this 
problem  by  searching  for  hyperbolic  trajectories  in  daily  averaged  fields.  Several  were  located  west 
and  south  of  Taiwan.  The  resulting  manifolds  were  then  superposed  on  the  hourly  EAS16  hourly 
velocities.  The  first  test  of  this  procedure  was  to  determine  if  particles  started  near  the  hyperbolic 
trajectories  followed  the  manifolds. 

Eigure  1  illustrates  the  findings.  It  shows  a  hyperbolic  trajectory  south  of  Taiwan  determined  from  the 
daily  averaged  EAS16  current  fields.  A  complete  animation  of  this  and  a  nearby  hyperbolic  trajectory 
is  available  at  http;//newark.cms.udel.edu/~mzweng/taiwan.  The  trajectory  and  infiowing/outflowing 
manifolds  were  located  in  these  fields  for  over  two  weeks.  In  Eigure  I,  particle  trajectories  and 
manifolds  are  superposed  on  hourly  model  velocities.  The  simulated  particles  shown  in  Eigure  1  were 
advected  by  the  hourly  fields  so  they  respond  to  small  time  and  space  scale  processes  that  are  filtered 
out  of  the  daily  averaged  fields  which  produced  the  manifolds.  Thus  this  calculation  represents  a 
reasonable  first  test  of  our  procedure.  The  movements  of  the  simulated  drifters  are  in  excellent 
agreement  with  the  manifolds.  We  did  not  observe  any  of  them  crossing  the  manifolds.  Moreover, 
they  generally  followed  the  evolution  of  the  manifolds.  The  next  step  is  to  compare  the  movement  of 
real  drifting  sensors  with  hyperbolic  trajectories  and  associated  manifolds  determined  from  the  daily 
averaged  velocity  fields. 

IMPACT/APPLICATIONS 

We  have  established  that  the  velocity  correlation  statistics  from  EAS16  generally  agree  with  similar 
statistics  from  deployed  drifting  sensors.  Surprisingly,  the  agreement  between  model  and  observed 
trajectories  is  poor.  One  obvious  reason  for  this  is  that  trajectories  of  real  drifting  sensors  are  very 
sensitive  to  subgrid  scale  processes,  which  are  only  parameterized  in  the  model.  The  remedy  for  this  is 
higher  model  resolution  and  assimilation  of  trajectory  data  in  EAS16.  However,  there  is  another 
possibility;  the  hydrodynamic  drag  of  drifting  sensors  causes  the  deviation.  Rectifying  this  issue  will 
require  detailed  knowledge  of  drag  characteristics.  Realistically,  one  expects  that  both  reasons 
contribute  to  the  model  trajectory  errors. 
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Figure  1:  Inflowing  (blue)  and  outflowing  (red)  manifolds  computed  from  daily  mean  velocities 
near  Taiwan  from  the  Navy  EAS16  model  for  the  period  2-8  April  2005.  Four  particle  positions, 
computed  from  hourly  model  velocities,  are  shown  as  colored  circles.  Hourly  EAS16  velocities  at  a 
depth  of  one  meter  (subsampled  at  every  fourth  grid  point)  are  also  shown. 
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